Colloidal gelation, percolation and structural arrest 
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By means of molecular dynamics, we study a model system for colloidal suspensions where the 
interaction is based on a competition between attraction and repulsion. At low temperatures the 
relaxation time first increases as a power law as a function of the volume fraction and then, due 
to the finite lifetime of the bonded structures, it deviates from this critical behavior. We show 
that colloidal gelation at low temperatures and low volume fractions is crucially related to the 
formation of spanning long living cluster, as we suggested in the last few years. Besides agreeing 
with experimental findings in different colloidal systems, our results shed new light on the different 
role played by the formation of long living bonds and the crowding of the particles in colloidal 
structural arrest. PACS: 82.70.Dd, 64.60.Ak, 82.70.Gg 



In attractive colloids a rich phenomenology is found in 
the temperature - volume fraction plane [l|-[3- At high 
temperature a hard sphere glass transition occurs at a 
volume fraction, <fi ~ 0.57. By decreasing the tempera- 
ture, at high volume fraction, the effect of the short range 
attraction produces an attractive glass line with a reen- 
trant line well described in the framework of the mode 
coupling theory 0, 0], and confirmed by experiments. 
At low temperature and increasing the volume fraction 
a much more complex situation arises H, char- 
acterized by a cluster phase followed by a kinetic arrest 
with viscoelastic properties very similar to those found 
in polymer gelation. 

One of the important problems is to understand the 
mechanism which gives rise to colloidal gelation. When 
phase separation does not occur, colloidal gelation has 
been interpreted in terms of cluster mode coupling the- 
ory or Wigner glass 0|. In contrast to these ap- 
proaches, we have suggested that at low volume fraction 
and low temperature, percolation phenomena play a cru- 
cial role in colloidal gelation characterized by the pres- 
ence of a spanning cluster of bonded particles [lj, [lj| 
as also found in recent experiments 18] . Of course, at 
high enough temperature and high volume fraction, one 
expects a structural arrest due essentially to the hard 
sphere part of the potential. At intermediate values of 
temperature and volume fractions, molecular dynamical 
studies have shown the occurrence of colloidal gelation 
phenomena characterized by a structural arrest, on aver- 
age well described by the mode coupling theory of attrac- 
tive glasses . Interestingly, the system is dynamically 
heterogeneous: The set of all particles can be partitioned 
in two subsets made of slow and fast particles and struc- 
tural arrest occurs close to the onset of a spanning cluster 
of only the subset of slow particles • One may ask the 
question whether the same scenario for colloidal gelation 
can be smoothly extended from this intermediate regime 



to the low volume fraction region. In this letter, by using 
a molecular dynamics simulations of a DLVO-type poten- 
tial for charged colloidal systems, we give evidence 
of a clear crossover from the low volume fraction regime 
to the intermediate regime. In the first regime we find 
the bond lifetime much larger than the structural relax- 
ation time. As a consequence, dynamically, the clusters 
behave as made of permanent bonds as in chemical gela- 
tion, and the regime is dominated by the dynamics of 
such clusters. Increasing </>, when these two time scales 
become comparable, the structural relaxation begins to 
be affected also by the crowding of the particles, and a 
clear crossover to a new glassy regime is found. Interest- 
ingly enough such crossover has also been found in some 
micellar systems at rather high volume fraction 0, , 
suggesting a common mechanism in both systems. 

We consider a system of N = 10000 particles, with 
(f> = 0.08-0.25, interacting via a DLVO-type potential, 
which contains a Van der Waals type interaction plus an 
effective repulsion due to the presence of charges: 
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where a\ = 2.3, a% = 6, 0,3 = 3.5, and A = 2.5 [2l|. 
With these parameters the repulsion term dominates the 
Van der Waals attraction at long range, providing a short 
range attraction and a long range repulsive barrier. The 
effective repulsion in the potential prevents the liquid-gas 
phase separation and stabilize the size of the clusters, as 
expected ^lj- The potential is truncated and shifted at 
a distance of 3.5a\ To mimic the colloidal dynamics, we 
performed molecular dynamics simulations at constant 
temperature. Equations of motion were solved in the 
canonical ensemble (with a Nose-Hoover thermostat) us- 
ing a velocity Verlet algorithm with a time step of 
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particles). We equilibrate the system at temperatures 
k B T = 0.2, 0.23, 0.25e and 4> increasing from 0.07 to 

0.23 m. 

Here we analyze the structure and the dynamics by 
varying the temperature and volume fraction. We calcu- 
late first the static structure factor, S(k), and the cluster 
size distribution, u(s), where two particles are linked if 
their relative distance is smaller than the local maximum 
of the potential [25| . 

At low temperatures, i.e. when the kinetic energy be- 
comes smaller than the repulsive barrier, and low vol- 
ume fractions, the static structure factor, S(k), displays 
a peak around fc — 2 (see Fig.^l. This feature is due to 
the fact that the competing short range attraction and 
long range repulsion produce stable clusters of a typical 
size. The presence of a typical size is also manifested in 
the cluster size distribution, n(s), plotted in Fig. [21 at 
ksT = 0.2e and 4> = 0.13. It is interesting to observe the 
secondary peak due to the fusion of two stable clusters. 
By increasing the volume fraction 4>, a spanning cluster 
appears at <fi = 0.16 for ksT = 0.2e, and the cluster size 
distribution typically displays a power law decay with 
an exponent t ~ 2.18, consistent with the random per- 
colation model [2^|. In Fig. [3] we show the dependence 
of the cluster size, s, on its radius, r, at 4> p : The data 
correspond to a compact structure for clusters of dimen- 
sion s < 10 and radius r < 1, a fractal dimensionality 
df ~ 1.2 on intermediate length scales, and a crossover 
to df ~ 2.5 at larger length scales. These results suggest 
that at low volume fraction compact stable clusters form 
with typical size s ~ 10 and radius r ~ 1. By increas- 
ing the volume fraction a residual attractive interaction 
between the clusters produces tube-like structures with 
fractal dimension df ~ 1.2 up a size s ~ 60. By further 
increasing the volume fraction these structures coalesce 
to build a random percolating network. For the temper- 
ature here studied the percolation threshold 4> p is weakly 
dependent on the temperature, 0.16 < <j) p < 0.17. 

A preliminary study of this model was published in 
Ref. [14J and these results were presented in the meet- 
ing "Dynamical arrest" , hold in Bad Gastein in Jan uary 
2005. This analysis has been lately developed in Ref. jlfij . 
where using a qualitatively similar interaction potential 
the authors have found again the behavior shown in Fig.s 
12 and |3J and stressed the role of percolation. Interest- 
ingly they have also shown that this model is actually 
able to reproduce the local order typically observed in 
the experimental system [l8|. 

In order to study the relaxation dynamics, we calculate 
different time autocorrelation functions at equilibrium. 
In Fig. 01 the incoherent scattering functions, F s (k,t), 
are plotted for k ~ 1.88, ksT = 0.25e and different <j>. 
From these functions we extract a characteristic relax- 
ation time, r a , defined by F s (k,T a ) = 0.1, plotted in 
Fig. El 

We also study the time correlation of the clusters in 




FIG. 1: The static structure factor, S(k), at ksT — 0.2e and 
4> = 0.13, displays a peak around feo — 2. 
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FIG. 2: Main frame: the cluster size distribution, n(s), 
at ksT = 0.2e and (f> = 0.13, displays a peak around an 
optimum cluster size. Inset: At <f> = 0.16 ~ <f> p , a power law 
tail appears. 



terms of the time autocorrelation functions of the bonds, 
B(t), which are plotted in Fig.0 
The function B(t) is defined as 



B(t) = 
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where riij(t) = 1 if particles i and j are linked at time 
t, Uijit) — otherwise. The data show that, differently 
from the behavior of F s (k,t), the bond relaxation does 
not change much below the percolation threshold, c/> p , 
and instead it becomes slower and slower above the per- 
colation threshold. By fitting the curves with stretched 
exponential functions, we obtain the bond lifetime 77,. In 
Fig. is compared with r Q . The data show that, at low 
volume fractions, the bond lifetime is much larger than 
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FIG. 3: The cluster size, s, as a function of the radius, r, at 
fcflT = 0.2e and <j> = 0.16 ~ <j> p : While at very small scales 
clusters are compact, on intermediate length scales the depen- 
dence suggests the presence of very thin and chain structures 
made up by clusters of typical size s < 10; It crosses to a ran- 
dom percolation type of structure over larger length scales. 




FIG. 4: (color online). The incoherent scattering function, 
F s {k,t), for k B T = 0.25e and = 0.08, 0.10, 0.12, 0.14, 0.15, 
0.16, 0.17, 0.19, 0.21, 0.23 (from left to right). Empty circles 
in figure correspond to values of <j> < <p p . 

the density relaxation time and in this region t q can be 
fitted by a power law. This result can be understood by 
considering that the system dynamically behaves as made 
of permanent clusters and the structural relaxation is di- 
rectly related to the size of the clusters, like in chemical 
gelation. At higher volume fraction close to the percola- 
tion threshold, the relaxation time and the bond lifetime 
eventually become of the same order of magnitude and 
one observes a deviation from the power law critical be- 
havior. 

Within this scenario, the percolation of long-living 
structures is the mechanism which produces the enhance- 



FIG. 5: (color online). The bond correlation function, B(t), 
for k B T = 0.25e and <f> = 0.08, 0.10, 0.12, 0.14, 0.15, 0.16, 
0.17, 0.19, 0.21, 0.23 (from left to right). Empty circles in 
figure correspond to values of <f) < cj> p . 
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FIG. 6: (color online). The bond lifetime, r&, and the relax- 
ation time, Tot, (in units of MD steps) as a function of the 
volume fraction cj> at k B T = 0.25e. The curve in figure is a 
power law fit, (cj> c — <f>) ~ 7 , with <f) c ~ 0.22 and 7 ~ 4.0. 

ment of the structural relaxation time at low volume frac- 
tion: Close to the percolation threshold, T a increases fast 
enough to reach the scale of the bond lifetime, and the 
dynamics crosses over to the glassy regime. 

In fact in this region on the time scale of the struc- 
tural relaxation time, the bonds cannot be considered 
as permanent anymore. This explains the crossover to 
a glassy regime, in which the structural relaxation time 
depends on the bond finite lifetime and the crowding of 
the particles. 

At higher temperature when the bond lifetime does not 
play any role, structural arrest is only due to the hard 
sphere component of the potential. 

By decreasing the temperature we find that the criti- 
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FIG. 7: (color online). The relaxation time, T a , as a func- 
tion of the volume fraction, (j>, for different temperatures, 
T= 0.25, 0.23, 0.2 (from right to left). The curves are power 
laws, {4>c — 4>) 1 > where respectively <f> c ~ 0.22 and 7 ~ 4.0, 
4> c ~ 0.21 and 7 ~ 4.3, and <j> c ~ 0.17 and 7 ~ 3.7. The per- 
colation threshold is weakly dependent on the temperature, 
0-16 < 4> P < 0.17. 

cal regime becomes more extended, due to the increase 
of the bond lifetime (see Fig. Q, and the critical volume 
fraction cf> c tends to the percolation threshold <j> p . There- 
fore at very low temperatures we expect the structural 
arrest to coincide with the onset of a spanning cluster on 
the time scale of the numerical simulations, in agreement 
with previous results and recent experiments^^]- 

In order to further understand and characterize the 
crossover region, we have calculated the dynamic suscep- 
tibility, defined as the fluctuations of F s (k,t) [2?| : 

X4 (M) = N((F s (k,t) 2 ) - (F s (k,t)) 2 ). (3) 

Similar quantities, largely studied in glassy systems (see 
Refs. |2ij|2j]), are related to the presence of dynamic het- 
erogeneities, spatial correlations between particles that 
move faster than the average. In Fig. Xi{k,t) is plot- 
ted as function of t, for k ~ 1.88, fc^T = 0.25e and dif- 
ferent values of <f). We find that Xi(k, t) exhibits a maxi- 
mum at time t* , which first increases and then decreases 
signalling the presence of a crossover as clearly shown 
in Fig. where the time, t* , and the maximum value, 
X&(k, t*), are plotted as a function of the volume fraction: 
At low volume fraction they both follow a power law di- 
verging, whereas a deviation from this critical behavior 
appears at volume fraction higher than the percolation 
threshold. 

The presence of a maximum at a time t*, increasing 
together with the height of the maximum was recently 
found in experiments on colloidal gels as colloidal gelation 
was approached [sfT ]. 

In conclusions in our model the presence of two well 
separated time scales at low temperature and low volume 



FIG. 8: (color online). The dynamic susceptibility, x*(k, t), 
for k ~ 1.88, k B T = 0.25e and <j> = 0.08, 0.13, 0.15, 0.17, 
0.21, 0.23 (from left to right). 
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FIG. 9: (color online). Main frame The maximum value 
of the dynamic susceptibility, x^(,k,t*), as a function of the 
volume fraction, cj>, at ksT — 0.25e. The curve in figure is 
a power law fit, (<f> c — 0) -7 , with <f> c ~ 0.17 and 7 ~ 0.68. 
Inset The time, t*, (in units of MD steps) as a function of 
the volume fraction, <j>, at fesT = 0.25e. The curve in figure 
is a power law fit, (<f> c — (f)) ' , with 4> c ~ 0.17 and 7 ~ 1.0. 



fraction originates a colloidal gelation related to the for- 
mation of persistent structures: Increasing the volume 
fraction the two time scales become comparable and a 
crossover from a gel-like behavior to a glass-like one is 
found. 
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